Spin conductivity in almost integrable spin chains 



o ■ 
O ' 

(N : 

<\ 
on : 



03 



i 

C 

O 

o 



> 
m 

00 

o 
o 



X 



Peter Jung and Actum Rosch 
Institute for Theoretical Physics, University of Cologne, 50937 Cologne, Germany. 

(Dated: February 1, 2008) 

The spin conductivity in the integrable spin-1/2 XXZ-chain is known to be infinite at finite 
temperatures T for anisotropies —1 < A < 1. Perturbations, which break integrability, e.g. a next- 
nearest neighbor coupling J', render the conductivity finite. We construct numerically a non-local 
conserved operator Jy which is responsible for the finite spin Drude weight of the integrable model 
and calculate its decay rate for small J'. This allows us to obtain a lower bound for the spin 
conductivity a s > c(T)/J' 2 , where c(T) is finite for J' — > 0. We discuss the implication of our result 
for the general question how non-local conservation laws affect transport properties. 

PACS numbers: 75.10.Pq, 02.30.1k, 75.40.Gb 



I. INTRODUCTION 

The behavior of transport properties of integrable 
systems have attracted considerable interest in the re- 
cent pasa^^^^^^^^^^il^^ . In such 
models due to the presence of conservation laws the cur- 
rents do not decay.— As a consequence, the dc conductiv- 
ity is infinite and characterized by a finite Drude weight 
D, Re ct(lu) — TtD5(uj) + a rcg (uj), where a lcg is the regular 
part of the conductivity. 

In real systems those conservation laws are violated 
by perturbations which often can be considered to be 
small. In these situations, the conductivity becomes 
finitei*2i2iii5iSi but remains very large as long as the per- 
turbations are small. 

A both theoretical 6 ! 7 ! 8 ! 9 ! 10 ! 11 ! 12 . 13 ! 14 , 21 and 
experimentall y 17 ! 18 ' 19 well studied example is the 
XXZ Heisenberg chain which is equivalent to a model of 
spinless Fermions with nearest neighbor interactions. 

The XXZ Heisenberg chain is integrable and therefore 
an infinite number of constants of motion exists for this 
model. All eigenstates can be uniquely labeled by a com- 
plete set of commuting operators, Q n with [Q n , Q m ] = 
for all n,m. The first two of these operators are given 
by the total magnetization Q\ = y~) k S% and the Hamil- 
tonian 



#0 = Q2=J2 h *' W 

i 

hi = J(SfS? +1 + SVSV +1 + AS?S* +1 ) (2) 

All other Q n can be constructed by a simple recursive 
formula^, Q n +i = [B,Q n ] with the so-called boost op- 
erator B = l/(2i) ^Zjjhj. All these conservation laws 
have a property which is important for the following dis- 
cussion: they are local operators in the sense that each 
Q n can be written in terms of a local "density" q n j at 
site i, 



E 



In 



(3) 



Besides the Q n , there exists a huge number of other- 
conservation laws Cj which can in principle be con- 
structed from the exact eigenstates \j), [H , = 
for Ej = Eji . In general, these operators are highly non- 
local objects in the sense that they cannot be written in 
the form for finite n. 

Only local conservation laws are associated with a con- 
tinuity equation d t q n ,k + 3n,k+i ~ jn,k = 0, where j n>k is 
the corresponding current density, and therefore only for 
local conservation laws a hydrodynamic description can 
be formulated. A main motivation for the present work is 
the question, to what extent non-local conservation laws 
are relevant in the sense that they lead to experimentally 
observable consequences in real materials. We therefore 
study the role of local and non-local conservation laws for 
transport in XXZ Heisenberg chains perturbed by weak 
next-nearest neighbor couplings J'. 

For J' — 0, both the heat conductivity n and the spin 
conductivity a s (or, equivalently, the electric conductiv- 
ity in the Fermionic language) are infinite and have a 
finite Drude weight. However, there is a main concep- 
tual difference between those two cases: the heat current 
Je is conserved (and actually given by Q$ as defined 
above), while the spin current J s does not commute with 
the Hamiltonian. Nevertheless, the presence of a finite 
Drude weight implies that a certain fraction of the spin 
current does not decay in time: part of the spin current is 
'protected' by conservation laws. This has been formal- 
ized many years ago by Mazur-22. and later generalized by 
Suzuki 23 . Suzuki showed that the Drude weight can be 
expressed in terms of correlators of the current with the 
conservation laws C;, 



(C?) 



(4) 



where q n ^ is local as it contains only spin operators 
on maximally n adjacent sites, i < j < % + n. 



where /3 = 1/T, N is the number of sites and the Ci have 
been chosen such that (CiCj) = for i ^ j. Note that in 
Eq. ([J) the sum runs over a basis of all conservation laws, 
local and non-local, commuting and non-commuting. 

Interestingly, it can be showni by simple symmetry ar- 
guments that the spin current is orthogonal to all known 
local conservation laws, {J s Q n } — for all n. Therefore 



2 



it seems that non-local conservation laws are responsible 
for the finite conductivity of the integrable model. What 
will happen to the spin current when the system is weakly 
perturbed, e.g. by a next-nearest neighbor coupling J' 
with J' <C J? For local conservation laws, e.g. the heat 
current Q%, the answer is known£ for small J', Q3 decays 
only slowly implying a large dc conductivity proportional 
to the life-time of Q3. Our present goal is to investigate 
whether the spin-conductivity shows a similar behavior. 

An example which shows that perturbation theory for 
local and non-local quantities can be drastically differ- 
ent has been discussed in Ref. 24|. In this paper it has 
been shown that an arbitrarily small inter-chain coupling 
can destroy a non-local order parameter (e.g. the string 
order of a spin-1 Haldane chain) in a gapped system. For- 
mally, the perturbations turn out to be proportional to 
the length of the system. In contrast, local order param- 
eters are always robust against small perturbations for 
all gapped systems. 

In principle one can try to investigate the transport 
properties for small J' directly by calculating the spin- 
conductivity from an exact diagonalization of the XXZ 
chain in the presence of finite J'. In such a calcula- 
tion Heidrich-Meisner et a/4 were able to show that the 
spin Drude weight vanishes in the thermodynamic limit, 
but a reliable determination of the resulting finite spin- 
conductivity is rather difficult even for large J'. Further- 
more, finite size effects grow rapidly^ for small J'. 

In the following, we will therefore use a different ap- 
proach based on a perturbation theory in J'. We con- 
struct numerically a non-local operator Jm which is con- 
served for J' = and responsible for the finite Drude 
weight of the unperturbed XXZ Hcisenberg chain. In a 
second step we derive a lower bound for the spin conduc- 
tivity of the perturbed system using results of Ref. [5(] and 
show that the spin conductivity is proportional to 1/J' 2 . 
Finally, we analyze to what extent Jm is a non-local oper- 
ator and discuss how the result can be interpreted. Ap- 
pendix [A] investigates the role of special values of the 
anisotropics where the Heisenberg model possesses extra 
symmetries. 



II. MODEL 

We consider the following Hamiltonian 

H = H + H U (5) 

where the XXZ Heisenberg chain Ho has been defined in 
Eq. (fTJ) and 

H 1 = J'Y J S%St +2 + StSt +2 + AS£Sf +2 (6) 
fe 

describes the (small) next-nearest neighbor coupling. For 
this model, the spin current J s is given by 



and we have omitted terms linear in J' as they give only 
subleading contributions to our final result. 

For J' = and — I < A < I, the Drude weight defined 
by 



Re a s (u)) = ttDS(lu) + <7 reg (w) 



(8) 



is finite^ ' 10 i 16 at T > as discussed above. Equation 
(HJ implies that the finite Drude weight is associated to 
constants of motion Ci of Ho with (CiJ s ) 7^ which we 
need to identify for our further analysis. More precisely, 
one can split the current operator into two pieces, 



with 



Js = J\\ + J.L, 



_ \ - {J s Cj) 

'ii-ElcfT 7 ''- 



(9) 
(10) 



Jy can be interpreted as the projection of the the spin 
current to the space of conserved quantities, i.e. the con- 
served part of J s and, indeed, one obtains directly from 
Eq. ©: 



D, 



(11) 



As described above, the known local conservation laws 
Q n do not contribute to J s , i.e. (J\\Q n ) = 0. J11 is a very 
complex non-local operator which is difficult to construct 
and handle analytically. For finite size systems with up to 
20 sites, however, one can construct J|| numerically using 
the exact eigenstates of Ho. As the Ci span the space 
of energy diagonal operators, we just keep the energy 
diagonal part of J,, i.e. 



= s 



E m E 



(12) 



For a finite value of the perturbation J' the Drude 
weight (|11[) is absent, as is known from numerical 
studies^iiS. which were, however, not able to investigate 
the regime of small J' due to large finite size effects in 
this limit. 

In Ref. [j| we have shown that a lower bound for the 
leading order contribution to the conductivity a a can be 
obtained in the limit of small J' by evaluating the corre- 
lation function T with respect to Hq: 

Rer(w) = - / «fte wt <[J||(t),J,|(0)]> . 



(13) 

As [Jm, H ] = 0, T(o-i) is proportional to J' 2 . The inequal- 
ity for the spin conductivity reads 



F(0) 



(14) 



where \ — P{J\\Js)/N = D s is the generalized (spin cur- 
rent) susceptibility and ^{ui)/x can be interpreted as a 
scattering rate of Jy, see Ref. [5| for details. Next we will 
present our analysis of the correlation function Eq. (Q2 
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FIG. 1: Leading order contribution to the spin current relax- 
ation rate for A = 0.75 and system size L = 20 for T — > oo. 
Finite size effects are small and T(lj) is finite at u) — as can 
be seen in more detail in the inset (thick line L — 20, dotted 
line L = 18, thin line L = 16, dashed line L = 14). 



FIG. 2: Scattering rate F(0)/x as a function of the anisotropy 
parameter A for L — 18, T = oo. For the isotropic system, 
A = 1, r(0) is zero, see text. The errors are comparable to 
the size of the symbols and are discussed in more detail in 
Appendix [X] 



III. NUMERICAL RESULTS 

We investigate T(u>) and the generalized susceptibil- 
ity x numerically in the T — ► oo limit via exact diago- 
nalization for system sizes up to L = 20 and for vari- 
ous anisotropics A using periodic boundary conditions. 
In this high temperature limit, the spin-spin correlation 
length vanishes and therefore finite size effects are small- 
est. Results for finite T> J (not shown) are essentially 
identical. 

The results for an intermediate A = 0.75 are shown in 
Fig.[TJ r(w) drops rapidly for small frequencies but satu- 
rates at a finite value. This saturation value lim^^o T(w) 
is almost independent of system size (see inset). This in- 
dicates that finite size effects are small despite the fact 
that J|| is expected to be a non-local operator. We there- 
fore conclude that for small J' 



TJ' 2 ' 



(15) 



in the thermodynamic limit. This is the main result 
of this paper: the spin-conductivity of a slightly per- 
turbed XXZ Heisenberg chain is very large, despite the 
fact, that the spin current is not protected by any local 
conservation law. For A = 0.75 we obtain for example 
c(T — > oo) = 0.92J 3 . For any finite temperature we ex- 
pect that the same result holds: in the limit of small J' 
the spin conductivity is proportional to 1/J' 2 . 

In Fig. [3 the behavior of the scattering rate T/x as 
a function of A is shown. Interestingly, the scatter- 
ing rate seems to vanish in the isotropic limit A — ► 1, 
r oc J' 2 (l — A) 2 . We have previously 6 observed the 
same effect for the scattering rate of the heat current, 
which turns out to be proportional to 1/J' 4 at the 
isotropic point. The reason for this unexpected result 



is that for the isotropic case one can construct an op- 
erator Q' 3 = Q 3 + J'AQ 3 such that the commutator 
[Q' 3 ,Hq + Hi] is of order J' 2 rather than linear in J'. 
As a consequence, the decay rate of the heat current at 
the isotropic point is proportional to J' 4 . Very likely, the 
same mechanism applies to Jn , too. A subtle and con- 
troversial issue d 10 ! 16 is the value of the Drude weight, 
D s = x, for A = 1. Both from numerics and from Bethe 
ansatz, there is evidence pointing either to a finite^ 6 - or 
vanishing 10,16 Drude weight in the thermodynamic limit. 
If the Drude weight vanishes for A = 1, our results are 
only of relevance for A < 1 . 

In appendix lAl we discuss a further effect: the Drude 
weight D s appears to be a discontinuous function of A as 
for special values of the anisotropies A = cos(7r/n), n — 
3,4,5..., one obtains different values for D s compared 
to anisotropies slightly away from these points. For the 
scattering rate T/x these effects are much smaller and 
possibly absent in the thermodynamic limit. 



IV. NON-LOCALITY OF J 

As stressed in the introduction, the spin current J s is 
orthogonal to all know local conservation laws Q n of the 
XXZ Heisenberg chain. This suggests that Jn, the con- 
served part of J s , is a non-local operator which cannot 
be written in the form of Eq. ([3]) . To quantify this state- 
ment, we expand the numerically constructed J|| in local 
operators A n i which contain products of spin-operators 
on n adjacent sites, 



<^> 1/2 ' 



(16) 
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FIG. 3: Relative weight Ysm=i Cm > see ^t G3> °f l° ca l op- 
erators with range up to n (n = 2, 4, 6, 8) contributing to 
as a function of inverse system size from N = 4 to N = 18. 
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shown in Figs. [3] and 0] describes which fraction of the 
operator Jm can be expressed in terms of operators of 
range n. For example, if one determines the c„ for H one 
obtains c 2 = J 2 /(J 2 + J' 2 ) and c 3 = J' 2 / ( J 2 + J' 2 ) . 
By construction one gets E n — p c n = 1 for a system with 
N sites. 

What types of behavior can be expected for c n ? First, 
one has to investigate whether c„ is finite or zero in the 
thermodynamic limit N — » oo. For example, for the 
square of a translationally invariant local operator (e.g. 
Hq), one finds that c„ drops proportionally to 1/-/V, such 
that limjv^oo c n — for all n > 0. Even if limjv^oo c„ 
is finite for each n, one can ask how rapidly lim^v^oo c n 
drops for n oo and whether E n =o um JV-»oo c„ equals 
1 or is smaller. 

As shown in Fig. [3l the c„ converge to finite values 
for N — > oo. For n = 2, this is a necessary consequence 
of the fact that the spin current is a range 2 operator 
and that the Drude weight of the spin current is finite. 
As the latter is proportional to (J S J||) 2 this implies that 
Jy has a finite overlap with a range 2 operator in the 
thermodynamic limit. 

A qualitative result of Fig. [3] is, however, that even 
operators up to range 8 have less than 40% of the total 
weight of Jy (but cs ~ 0.02 is already very small). As 
E^Li c-n, — 1, the c„ have to drop faster than 1/n for 
large n in the thermodynamic limit. Fig. |4] shows that 
the c„ decay extremely slowly with n. In this sense Jy 
appears to be a rather nonlocal operator but we cannot 
decide from our numerics whether E^Lo um A r -+oo c„ = 1 
or smaller. 



range 



FIG. 4: Relative weight c n of local operators of range n (n = 
3, ... 8) contributing to Jm for different system sizes. Note that 
in the thermodynamic limit most of the weight is carried by 
operators involving more than 8 consecutive sites (see Fig. [3]). 



where the A n i define a complete orthogonal basis in 
the space of operators, (A n iA m j) = for n ^ m or 
i =/= j. The A n i are written as sums of products of spin- 
operators, where each product contains spins on n adja- 
cent sites. Here we use - as above - the (T = oo) expecta- 
tion value as the scalar product in the space of operators. 
In Eq. (|16p obviously only translationally invariant her- 
mitian operators contribute which also conserve S z . For 
n = 1 there is just one such operator, An — J2i$ii 
for n — 2 one finds 3 such terms A 2 \ = J2i SfSf +1 , 

A 22 = Y,iSt S i+i + h - c -i A 23 = iEj^S^i - h.c). 
The 10 operators of range 3, A 3i , contain both products 
of two spin operators, e.g. S*Sf +2 and products of 
three spin-operators, e.g. S* Sf +1 S* +2 - 



V. CONCLUSION 

In this paper, we have shown that the spin- 
conductivity of a one-dimensional anisotropic spin-chain 
is strongly enhanced close to the integrable point. It di- 
verges (at least) as 1/ J' 2 for J 1 — > 0. This is the expected 
behavior for a situation where a local conservation law 
prohibits the decay of the current at the integrable point. 
However, as emphasized by Zotos, Naef and Prelovsek-, 
the spin-current is orthogonal to all known local conser- 
vation laws of the XXZ chain. 

There are two possible interpretations of this result. 
First, the conserved part J\\ of the spin current could 
nevertheless be 'sufficiently' local to define a slow hydro- 
dynamic mode. Second, the theoretical prejudice, that 
only local conservation laws (i.e. those associated with a 
continuity equation) lead to slow modes, may be wrong. 
In this respect, the results of section|lVl where this ques- 
tion is investigated, are ambiguous. On the one hand, 
we could prove that J|| is a highly non-local operator in- 
volving products of operators acting on widely separated 
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FIG. 5: Drude weight D = \ (upper panel) and scattering 
rate r(0)/x (lower panel) for the anisotropy A = cos(7r/3) = 
1/2 and a nearby value A = 0.499 as a function of inverse 
system size 1/N. The error bars represent the uncertainty 
inherent in the fitting procedure. 
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Forschungsgemeinschaft through SFB 608. 



APPENDIX A: SPIN-CONDUCTIVITY CLOSE 

TO AND AT A = cos tt /n 



In this appendix we discuss the behavior of the spin 
conductivity for anisotropics A = cos(n/v). At these 
special points it is known that there are further sym- 
metries which, for example, simplify the Bethe ansatz 
equations considerabl y 25 ' 26 . Interestingly, at these spe- 
cial points thermodynamic quantities show unexpected 
logarithmic corrections^ 6 .. 



sites. On the other hand, the relative weight of range-n 
operators, c„, is finite in the thermodynamic limit. 

In this paper we have shown that the transport proper- 
ties of simple one-dimensional problems depend quanti- 
tatively and qualitatively on 'exotic' and rather complex 
conserved quantities. For the future, it would be inter- 
esting to gain a more analytic understanding of these 
conservation laws. 
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In Ref. 0], Naef and Zotos found numerically that 
the Drude weight at these special points differs for fi- 
nite systems significantly from the values obtained for 
slightly different anisotropies. They concluded, how- 
ever, that these differences vanish in the thermodynamic 
limit. While we have reproduced the numerical results of 
Naef and Zotos, we interpret our results differently. In 
Fig. [5] we show the Drude weight of the integrable model, 
D s = x, and the scattering rate, T(0)/x , as a function 
of 1/N both for A = costt/3 = 1/2 and A = 0.499. Ap- 
parently, different values are obtained for N — > oo for \ 
(subtle logarithmic finite size effects can possibly invali- 
date this analysis) while the effect for the scattering rate 
r(0)/x is much smaller (and possibly absent). 
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